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Abstract 

We study the magnon modes in the presence of a topological soliton in a 
2d Heisenberg easy-axis ferromagnet. The problem of magnon scattering 
on the soliton with arbitrary relation between the soliton radius R and the 
"magnetic length" Ao is investigated for partial modes with different values 
of the azimuthal quantum numbers m. Truly local modes are shown to be 
present for all values of m, when the soliton radius is enough large. The 
eigenfrequencies of such internal modes are calculated analytically on limiting 
case of a large soliton radius and numerically for arbitrary soliton radius. It is 
demonstrated that the model of an isotropic magnet, which admits an exact 
analytical investigation, is not adequate even for the limit of small radius 
solitons, R <C Ao: there exists a local mode with nonzero frequency. We use 
the data about local modes to derive the effective equation of soliton motion; 
this equation has the usual Newtonian form in contrast to the case of the 
easy-plane ferromagnet. The effective mass of the soliton is found. 

75.10.Hk, 75.30.Ds, 75.70.Kw 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Nonlinear topologically nontrivial excitations (solitons) are well-known to play a spe- 
cial role in low-dimensional magnetic systems. For example, the presence of vortices in 2d 
easy-plane (EP) magnets gives rise to the Berezinskii-Kosterlitz-Thouless (BKT) phase 
transition.^!! Kinks in Id magnets and localized Belavin-Polyakov (BP) solitonsi in 2d 
isotropic magnets are responsible for the destruction of long-range order at finite tempera- 
ture. The soliton signatures in dynamical response functions can be observed experimentally. 
Translational motion of solitons leads to the so-called soliton central peak. Another pos- 
sibility to detect soliton signature is to look for magnon modes, localized at the soliton (local 
modes, LM); such modes are the cause of soliton magnetic resonance at the characteristic 
frequencies of "intrinsic" motion.0 

This can be explained within the scope of the so-called soliton phenomenology, where 
the magnet can be described as a two-component gas of elementary excitations: solitons 
and magnons. Such an approach was developed for Id magnets.! It was shown that the 
contribution of magnons and soliton-magnon interaction is important for this approach, 
as is obvious, for example, in the discussion of soliton magnetic resonance. I But for Id 
magnets the soliton-magnon scattering causes a change of the magnon density of states, 
which is necessary for a self-consistent calculation of the temperature dependence of the 
soliton density.0 

For the 2d case, the concept of a soliton (vortex) gas has been extended to describe 
EP magnets above the BKT vortex-antivortex unbinding transition, with a finite density of 
vorticesJE The concept of localized topological solitons (^-topological solitons, see details 
in Ref. [[!]) was used to explain the EPR line-width in easy-axis (EA) magnets.EUll The 
classical example of a localized soliton is BP-solution, which exists in isotropic magnets. 
The scattering problem for such solitons allows an exact analytical solution because of the 
scale- and conformal- invariance of the model.[iiB For EA magnets it is necessary to take 
account of some important modifications due to the breaking of such symmetries. That is 
why, both for EA and EP 2d magnets, the soliton density has not been calculated, but has 
been used as an input parameter. The general features of the 2d soliton dynamics are not 
clear at present. In particular, the form of inertial terms in the dynamical equations for the 
soliton centre is unknown; so the study of localized states is an issue of current research. 
Note that the problem of the soliton dynamics and the problem of the existence of LM 
are intimately connected to each other, and to the problem of soliton-magnon scattering. 
For example, using numerical data for the scattering amplitude, a non-Newtonian effective 
equation of motion of the magnetic vortex was constructed.^ 

Usually non-Id solitons and, especially, the problem of soliton-magnon scattering, are 
treated numerically.0HHl Only in the case of purely isotropic magnets, with the exact BP- 
solution, the analytical investigation of the scattering problem has been done recently.Bi!§ 
In particular, it was found that a soliton has a number of local modes of zero frequency;il 
such modes determine the main features of the scattering pictured However, even a small 
uniaxial anisotropy leads to principal changes. Truly local modes (in contrast to quasi-local 
ones, below) with nonzero frequencies were found numerically in Ref. |26] for the soliton in 
EA ferromagnets (FM) for some special values of the azimuthal quantum number m; namely 
the mode with m = — 1 for the soliton with the topological charge q = 1, and the mode with 
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m = 2 for q = 2 were investigated. Note that local modes do not exist for the models of 
easy-plane and isotropic ferromagnets, where the magnon spectrum has a gapless dispersion 
law. The case of the 2d EP antiferromagnet is an exception; for this model a finite-frequency, 
truly localized internal mode with m = exists inside the continuum spectrum.0 

Recently there has also been renewed attention to the internal dynamics of the soli- 
ton. It is caused by the studying of small ferromagnetic particles in the so-called vortex 
state (magnetic dots). Magnetic dots are micron-sized magnetic samples, placed on a non- 
magnetic substrate; they have different shapes: circular, quadratic, etc. The magnetic dots 
are interesting for the usage in the high-density magnetic storage devices@i! and from a 
fundamental point of view. As the dots are rather large (the size of magnetic dot exceeds 
the critical size of a single domain), their ground state could be inhomogeneous. Such a 
dot in the vortex state contains a magnetic singularity in the centre of the vortex, or the 
Belavin-Polyakov soliton type.0 The discrete magnon modes were observed for the uni- 
formly magnetized dots by a resonance technique;^ their theory was constructed using the 
Landau-Lifshitz equations!!. 

In this paper we consider the magnon modes which exist in 2d Heisenberg EA FM with 
a topological soliton. The model of an easy-axis ferromagnet is presented in Sec. ||. An 
approximate solution for the soliton structure is examined. In Sec. |TTT] we formulate equa- 
tions, describing magnons and solitons. In Sec. truly local magnon modes are calculated 
analytically by variational approach for the limiting cases of large and small radius solitons, 
and numerically for arbitrary radii. In Sec. |V| the problem of soliton-magnon scattering 
is formulated and treated numerically for different partial modes. The phenomenology of 
a pseudopotential is proposed in Sec. |V| to analyze analytically both local states and the 
scattering problem. In Sec. [VI I| we use results for local modes to derive the equations of 
soliton motion. Such equations are of Newtonian type (with the account of gyroscopic force), 
in contrast to both the case of the BP-soliton in the isotropic magnetic and the case of the 
vortex in the EP mag In the framework of this equation the soliton effective mass is 

calculated. A discussion and concluding remarks are presented in Sec. [VIII. 



II. THE MODEL AND ELEMENTARY EXCITATIONS 

The dynamics of the classical ferromagnet is described by the Landau-Lifshitz equations 
for the normalized magnetization m. In angular variables, m x + im y = sin9 exp(i<ft), these 
equations correspond to the Lagrangian 

L = Aj rf»a; |I(i-oos0)|£-W(M)} (1) 
with the energy density A ■ W(9, <fi), 

<j>) = \ ( V9f + 1 ( V0) 2 sin 2 9 + ^ sin2 ( 2 ) 

where A = JS 2 , J is the exchange integral, S the atomic spin, Ao = \J A/K the charac- 
teristic scale ("magnetic length"), K the energy of the anisotropy, and D the spin-wave 
stiffness. 
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The model of EA FM has well-known magnon excitations 9 = const < 1, <f> = ut — k-r 
above the homogeneous ground state. The dispersion law for the magnons has a finite 
activation frequency u> = D/Aq, 

u(k) =u + D\k\ 2 , (3) 

where k is the wave vector. 

The simplest nonlinear excitation in the system is the dynamical (precessional) topolog- 
ical solitonlHl 



<t> = <fo + QX + = 6 (r), 6 (0) = n, 6 (oc) = 0, (4) 

where r and \ are polar coordinates in the plane of the magnet, the integer q plays the role 
of a topological charge, and Q (0 < Q < Uq) is the frequency of the internal precession. The 
necessity to consider solitons with internal precession is caused by the fact that, according 
to the Derrick-Hobart theorem, static solitons in models like are unstable. Formally, 
there is no function 6 (r), which could provide a minimum of the energy 





E = nAl rdr{[^) +stfe [^+?)}. (5) 

The precessional soliton corresponds to the conditional minimum of the energy (Bp for a 
given number N of magnons bound in the soliton 



N, 



AT=_1/ rdr (1 - cos o (r)) , (6) 
Jo 

which is conserved for the uniaxial case (A^2 = 27rA/hujQ is the characteristic number of 
bound magnons in the 2d soliton) ,0 see also Refs. PIJ53L This corresponds to the condition 
5 [E — NhQ] = 0. The form of the function 8o(r) is defined by the ordinary differential 
equation, 

d% ld6 ( 1 q 2 

H ; sin V cos (J Q [ + — 



dr 2 r dr \Aq r 2 

u A 2 



+ ——2 sin O = O, (7) 



with the boundary conditions The precession frequency Q is fully determined by the 
number of bound magnons N. The dependence Q(N) was calculated numerically in Ref. |34] ; 
large radius solitons (N 3> N%) have a very low precession frequency (Q ujq); and Q — > ujq 
for small radius solitons (N N2). 

In the case of large radius solitons, R ^> Ao, the approximate "domain wall" solutionS 

r - R 

cos0o( r ) = tanh — - — (8) 
Ao 

is applicable, which describes a Id domain wall at the distance R from the origin. One can 
obtain the relation between the soliton radius R and the frequency of the soliton precession 
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£7 from variational condition, in which R plays the role of a fitting parameter for the function 
(§). In the first approximation Q(R) m ujq ■ A /R. 

The opposite case of small radius solitons, R -C A , corresponds to the Belavin-Polyakov 
solution tan(6 l /2) = (R/r)^, which describes the isotropic limit of the model ([D; it is valid 
for distances r <C A J At larger distances the algebraic decay is replaced by the general 
exponential form, 



O (0 oc exp (-r/r ) , r = W — — — . (9) 

The dependence Q(R) was calculated both analytically and numerically in Refs. |34| , |35| . The 
frequency of the soliton precession satisfies Q — > ujq when R — > 0, but the dependence f2(iV) 
is not analytic, and dVt/dN — > oo as N/N 2 — ^ 0.0 

To describe approximately the soliton structure in the intermediate case of arbitrary R, 
we use a fitting method.il We choose a trial function of the form 



0(fit)/ r) ^(fit) 
tan — - — = • exp 



(r - R^) 

r 



(10) 



which has valid asymptotics for r — > and r —>■ oo (we discuss the case q — 1 only). The 
parameter i?( fit ) can be found numerically by fitting the trial function (fTOj ) to the numerical 
solution of Eq. (0). Such an approximation shows that the trial function describes the soliton 
shape with an accuracy of about 10~ 2 , see Table Q. From this table one can see that the 
values of the fitting parameter R( &t > are very close to the "exact" soliton radii; it should 
be understood by the fact that 6q (r = R^j — ir/2. Numerical data for 9$(r) and fitting 

results for 9q (r) are plotted in Fig. [I]. 



III. MAGNON MODES 

To analyze magnons on the soliton background, it is convenient to introduce local coor- 
dinates {ex, e 2 , e 3 }, which describe the configuration of the magnetization unit vector m in 
the unperturbed soliton: e 3 coincides with m in the soliton (f|), e.\ = e^cos^o — e^sin^o, 
and e.2 = e% x e.\. Then the linear oscillations around the soliton can be described in terms 
of projections of the magnetization m on the local axes: /i = m ■ e\ and $ = —m • e 2 . In 
the absence of the soliton, # — 0, 0o — Qt, such oscillations correspond to free magnons in 
the rotating coordinate frame. 

The linearized equations for $ and \i can be presented in the form: 

r 2 . ,n 2qcos6 dfi 1 dfi 

[" v +Ul ^ + ^^d^ = ^W (lla) 

r „ 9 , .-, 2(7 cos On d'd 1 dfi x 

I- v +u ^"-^i* = -^m (ub) 

with the "potentials" 
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U x {r) 
U 2 (r) 



A 2 ' r 2 



n 



+ ^ I cos 26> (r) ^ cos 9 (r), 



de (r) 



(12a) 



dr 
d9 (r) 



dr 



+ cot 6 {r)V% 



+ ( ^ + ^ ) cos 2 W 







cos (r), 



(12b) 



where V 2 is the radial part of the Laplace operator. 

It is convenient to represent the solutions of Eqs. ([11]) in terms of a partial-wave expansion 
for t? and /i 



E 



u n (r 



v a (r)) sm(mx + uj a t), 



(13a) 
(13b) 



where a = (k, m) is a full set of eigennumbers, k and m being radial and azimuthal quantum 
numbers, respectively. The quantity uj = uj — Q is the magnon frequency in the rotating 
frame; k and u are connected by the dispersion law, cf. Eq. (||), 



ap 



D\k\ 



UJ, 



gap 



(i) 



where u> gap is the magnon gap frequency in the rotating frame. Below we will replace uj by uj, 
almost everywhere without confusion, because we will work in the rotating frame mainly; the 
index a will be omitted, too. Thus we then finally obtain the following eigenvalue problem 
(EVP) for the functions u and v: 



H x u = 
H 2 u = 



m „/ \ ^ 
+ -r + V(r) - 



uj AI 



1171 t r / \ UJ 

H + — - V(r) + 



UJqAI 



u = W(r)v, 
v = W(r)u, 



(14a) 
(14b) 



where Hq = —V 2 . + Uo(r) is the 2d radial Schrodinger operator; the "potentials" are 



U n (r) 



(Wi+W 2 



V(r) 



2qm cos 9 



(15a) 



Note, that the Eqs. (0) are invariant under the conjugations uj — ► — uj, m —>■ —m, and 
u <-> v. In a classical theory we can choose any sign of the frequency; but in order to make 
contact with quantum theory, with a positive frequency and energy Sk = h\uJk\, we will 
discuss the case uj > only. Thus, there are two different sets of Eqs. (|T4] ) for m — |n| and 
m = —\n\. Besides there appears a difference between Eqs. (14a) and (|14b|) . Since ( 14b ) 
has the asymptotically equivalent form, — V^f + ujv = 0, the function v has an exponential 
behaviour; at the same time Eq. (|14a| ) yields oscillating solutions. So, by choosing uj > 0, 
the variable v becomes a slave variable in the EVP (|l"4T). 
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The "coupling potential" W, 



W{r)= l -{U 2 -U l ) 



sm 0q[ a| + ^J-^ 



(15b) 



is positive, it has a maximum value at a distance of about the soliton radius R from the origin, 
vanishing both at r = and far from the soliton. Using the "domain wall" approximation 
(§), one obtains that the maximum value of W oc 1/-R 2 , when R ^> A . In the opposite case 
of small R, in the range r < R the soliton looks like the Belavin-Polvakov solution; here 
r(d9o/dr) « — sin#o ; an d W = in the main exchange approximation.!^ Therefore, one can 
suppose that W <^1/Aq for any R. In fact, this is confirmed by numerical calculations, see 

Fig. @. 



However, one can see below in Sec. [TV] that the problem is impoverished when the "cou- 
pling potential" W(r) is ignored; even in the limits of small (large) soliton radii, the main 
properties of local modes vanish, if we assume that W(r) = 0. The reason is that a set of 
more than one coupled Schrodinger-like equations has properties, which are absent from a 
single equation. To explain this fact, remember that for a single Id Schrodinger equation 
there exist two independent solutions; for the eigenfrequencies inside the continuous spec- 
trum such solutions have oscillating asymptotics. For a set of 2, 3, . . . equations there exist 
4, 6, ... linearly independent solutions, respectively; some of them can have exponential 
asymptotics. Naturally, the real modes have an oscillational form here; we will use this fact 
below for the numerical analysis. But it is necessary to take into account the exponential 



solutions, too. In particular, as was shown in Ref. |27|, such exponential asymptotics could 
correspond to truly localized states inside the continuous spectrum, which are forbidden for 
equations of the Schrodinger form. In the reduced EVP with two independent equations 
like ( |i~4"D with W = 0, there is a continuous spectrum uj > uj gap for the function u(r) (u ^ 0, 



v = 0) and u < — uj gap for the function v(r) (u = 0, v ^ 0). There are no limitations for 
the lowest value of the discrete levels; in particular, the frequencies of the discrete spectrum 
could be less then — a; gap for the function u. At the same time, as follows from the general 
set of coupled Eqs. (0), frequencies of local modes lie in the region — u gaiP < uj < u; gap . 
Therefore, it is necessary to solve the general EVP fli~4]) for an exact analysis of the problem. 

We should note here that the model of an isotropic magnet is a special case, in which 
the "coupling" potential is exactly equal to zero. To study an EVP like ([14]) with W — 0, 
it is sufficient to take into account one equation only (namely, the equation for the variable 
u); the second one, which has unphysical solutions for uj > 0, was ignored in Ref. [19|. That 
is why the results for EA FM, obtained in this article, and the results for the isotropic 
magnet^ differ strongly even for small anisotropy. One difference is caused by the presence 
of truly local modes. Another one is due to the peculiar structure of the quasi-continuous 
spectrum with uj > u; gap for a finite geometry: the spectrum of the EA FM has a system of 
doublets with opposite signs of m, while there are no doublets for the isotropic magnet. 
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IV. LOCAL MODES 



According to the dispersion law (|3[), free magnons only exist when uj > uj gap . It is 
naturally to look for local modes in the range < uj < co> gap (remember that we discuss the 
case D > ). In principle, such local modes could exist inside the continuous spectrum, as 
for the case of the nonlinear cr-model, describing the EP antiferromagnet.0 However, our 
study shows that this is not the case for the EA FM. Existence of local modes in such a 



magnet was found numerically in Ref . ^6] for a soliton with some special parameters; namely 
the mode with m = — 1 for the soliton with the topological charge q = 1, and the mode 
with m = 2 for q = 2 were calculated. Such modes appear in field theories, which include 
soliton-like solutions.!! 

To study local modes, note that near the soliton centre, r <C Aq, one can obtain the 
asymptotics of the magnon amplitudes in the following form: 

(\ \i~ m \ 
— ) [l + O(r 2 /A 2 )}, (16a) 

(\ |<j+rrt 
[l + O(r 2 /A 2 )}. (16b) 

Here the coefficient A can be determined by the normalization; the presence of a nonzero 
factor e is caused by the "interaction" W(r) between Eqs. ( |14a| ) and ( |14b| ); its value cannot 
be found through this asymptotic expansion. 



In the opposite case of large distance, r 3> Ao, the EVP (14) has the following gen- 
eral solution (here and below we consider the most interesting case of solitons with unit 
topological charge q = 1, which has lowest energy) 

u m (r) 7H A m K\ m+1 \(>c^r) + B m I\ m+1 \{jc[r) (17a) 
v m (r) w C m K\ m _ x \{x1[r) + D m J| m _i|(x^r), (17b) 

where K v and I u are Macdonald and modified Bessel functions respectively, and xf = 
\J (u>o — Vl±uj) / D. Note, that at large distances, z 3> 1, the modified Bessel functions 
have exponential behaviour, I u (z) oc exp(z)/-y/i, K u (z) oc exp(— z)/y/z, so we use the 
asymptotical condition u m , v m oc exp(— xfr)/y/r to obtain the localized solutions. 

Firstly note that for the special cases m = 0, +1 there exist "zero-modes" with u = 0. 
Such modes are caused by the internal symmetry of the problem. They have the same form 
as for EP magnets. One of the zero-modes, the so-called translational mode with m = +1, 

d0 sin 6» d6 sm6 

u+x = —, , v +1 = — H (18a) 

dr r dr r 

describes the position shift of the soliton. Another one with m = 

u = sin# ) t>o = — sin^o (18b) 

is a rotational mode, which corresponds to the presence of an arbitrary parameter, the phase 
ipo, in the soliton structure (f|). For the isotropic magnet, as well as for the EP magnet, 
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similar modes exist, too; but there they are quasi-local due to their slow power-law decay, 
1/r, or 1/r 2 . In the EA magnets the decay has an exponential form. 

The remaining modes could be investigated numerically; analytical results were ob- 
tained for some limiting cases only, see below. To study the EVP numerically, we use 
the two-parameter shooting method for the integration of the set of Eqs. flU]); a similar 
one-parameter shooting scheme was proposed in Ref. |27] to study the vortex-magnon scat- 
tering problem in 2d antiferromagnets. Numerically, we "kill" two growing exponents in 
(|17|) by choosing two shooting parameters e and uj. As a result, we have found numerically 
that there are a number of local modes. Modes with opposite signs of m bind to doublets; 
this picture is well-pronounced for \m\ > 1 in the limit of large R, see Fig. |3]. 

Such local modes exist only for not too small soliton radii (not for large Q). Eigenfre- 
quencies for such modes with arbitrary azimuthal numbers m (|m| > 1) increase rapidly with 
the decreasing of the soliton radius, and reach the boundary of the continuous spectrum. It 
looks as if all such modes "leave" the region of the gap in the spectrum, if the soliton radius 
is smaller than some critical value R c m , see Table ||. It is natural to suppose that all modes 
with \m\ > 1 transform to the quasi-local ones with decreasing of R; this fact can be tested 
numerically due to the singularities of the scattering amplitudes on such quasi-local modes, 
see the next section. 

Thus, in the range R < 1.52A there is only one local mode with non-zero frequency, 
m = —1, see Fig. |3]. That is why only the mode with m = — 1 was observed in Ref. |26| for 
the soliton with q = 1. 

The eigenfunctions of bound states u and v are localized near the soliton core. The 
functions u and v are plotted in Fig. f| for the mode with m = — 1; the picture of well- 
localized eigenfunctions is specially pronounced for the large radius solitons, see Fig. ||(a). 

In order to analyze the spectrum of local modes analytically, let us reformulate the EVP 
(|14]) as a variational problem for the functional 

F[u, v] = J rdr (uHiu + vH 2 v - 2Wuv ) . (19) 

We start with large radius solitons, R ^> Ao, when the "domain wall" approximation (§) 
is valid. In the region of interest, r ~ R ^> A , the operator H from the EVP ( |T^D has in 
the main approximation for A /R a simple reflectionless potential, 



cosh 2 ((r-i?)/A ). 
with a single localized state tpo, corresponding to the zeroth eigenvalue 

U = 0, Mr) = (2RA y 1/2 I cosh ((r - R)/A Q ) . 

The other potentials in the EVP can be considered as corrections to L; this is the reason 
to put the trial functions in the form u = aipo, v = bipo with trial parameters a and b. 
The variational approach with such simple trial functions leads to eigenfrequencies for the 
different modes in the form 

u m = ma + ^/(m 2 /?-<5) 2 - 7 2 , (20) 



dr 2 + A| 
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V'o 



2 cos 6n(r) 



l=(^\W(r)\^ ), 5=(^ 



1 d 

r dr 



V'o 



+ Ocos 9o(r) 



-0C 



where (ipo \. . .\ i/jq) = J °° ipo(. . .)ipordr . However, there is a difference between the real 
solutions u, v, and the trial function ip . This difference provides, e.g., the existence of 
"zero-modes". We do not like to improve the form of the eigenf unctions, but only the 
coefficients in the Eqs. (20), which describe the eigenfrequencies. Using the additional 
condition for "zero-modes", 



L) r , 







for m = 0, +1, 



(21) 



the eigenfrequencies uj m can be expressed through two independent coefficients only. It is 
convenient to present the eigenfrequencies in the form 



uj m = ma + \m\\/a 2 + (3 2 (m 2 — 1). 



(22) 



Using trial functions ipo, we can calculate the parameters a and j3 for the large radius 
solitons, 



a 



-u (A /R) 3 
In this case the spectrum has the form 



P = u (A /R) 2 . 



^0 




Ao 
R 



(23a) 



(23b) 



This formula corresponds very well to the numerical data for R 3> A , see Fig. ||. In the 
main approximation for Aq/R there is a system of doublets (for all modes with |m| > 1) 
with quadratic dependence of the mean frequency on A /R, 



(Jo\m\\/m 2 



R 



(23c) 



and the splitting is a small cubic correction, Auj± m oc m(Ao/R) 3 . For the special case 
m — — 1 the dependence u(R) has the form 



U) m= -i 



2Un 



Ao 
R 



(23d) 



Eq. ( P3b|) describes qualitatively all modes with arbitrary R. In particular, it describes 
the fact, mentioned above, that modes reach the boundary of the continuous spectrum at 
some finite frequency of the soliton precession Q^. Let us estimate the positions of the 
crossover points, where the modes "leave" the discrete spectrum. Using Eq. (|23b| ) , one can 
obtain the condition for the critical values of x = VL c m /ujQ\ 



\m\ a/x 6 + (m 2 — l)x A — mx c 



1 — x. 



(24) 



10 



Results are very close to the numerical data, see Table [T| We expect that the local modes 
transform into the quasi-local ones inside the continuous spectrum; but their frequencies 
can still be described by the formula ( [23 b| ), see Sec. [VT| . 

For the small radius solitons, R -C A , there only exist local modes with \m\ = 1. Thus, 
we will consider this case. Constructing the trial functions, one can suppose that they should 
be similar to the modes with m — ±1. Eigenfunctions for the translational "zero-mode", 
u + \ and f+i, are well-known, see Eq. ( |18a| ). To construct the trial functions it is convenient 



to use functions w_i and v _i in the form 

u_! = r 2 u +1 , v_i = r 2 v +1 , 

which coincide with the exact "zero-mode" eigenfunctions for the mode with m = — 1 in the 
limiting case of the BP-soliton.El Now we can choose the trial functions u and v as follows, 

u = a ■ u + i + b ■ v = c ■ v + i + d ■ V-x. (25) 

The variational problem ([19]) can be solved analytically in the main approximation for R/A 
using the explicit BP-solutions. Tedious calculations lead to the eigenfrequency for small 
R, 

u m= -x « C ■ {u - Q) . (26) 

Unfortunately, this method is not exact and leads to a value of the constant C which is 
larger than the numerical value. This is caused by the fact that the trial functions, which are 
constructed from zero modes, have the wrong asymptotical behaviour (Q) at large distance, 
6* (r), 0' (r) oc exp(— r/r ), instead of (|17|), i.e. u oc exp(— xf r), v oc exp(— x^r). To improve 
this result, one can change the asymptotical behaviour of the functions ( |23| ) by exponential 
factors: 

u — >u ■ e r (^ r °-^), v — >v ■ e r ( 1/r °-^). 

In this case the calculations can be done numerically only; the value of the constant C ~ 1. 
For the approximate description of the eigenfrequency of the mode with m = — 1 in the 
intermediate case of arbitrary R, we use the trial function 

(trial) (n _ 2 C 3 (1 ~ C) , _ / 9? x 



Providing the asymptotically correct behaviour, ( |23d[ ) and (PB|), this approximation repro- 



duces to the numerical results with an accuracy of about 2-10 3 , see Fig. [| 

V. THE SCATTERING PROBLEM: NUMERICAL RESULTS 

Let us describe the scattering of magnons by a soliton. Note that without a soliton, free 
magnon modes have the form 

^m(^) oc J m (kr ), v m (r) = 0, (28) 
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where A; is a "radial wave number", k and u are connected by the dispersion law (^'), and 
J m are Bessel functions. The free-modes like u m play the role of the partial cylinder waves 
of a plane spin wave 

oo 

exp (ik-r- tout) = i m J m {kr)e imx - iut . (29) 

m=— oo 

To describe magnon solutions in the presence of a soliton, one should note that far from 
the soliton, r>Ao, the "coupling potential" W(r) is exponentially small. Therefore there 
are two asymptotically independent EVP; the magnon amplitude v m has the same form 
(ELZlD as f° r ^ e l° ca l modes, but the function u m shows an oscillating behaviour like 

u m (r) oc J| m+ i|(b~) + a m (k)Y\ m+1 \(kr), (30) 

where Km+il are Neumann functions. The quantity cx(fc) derives from the soliton-magnon 
scattering; it can be interpreted as the scattering amplitude. 

We use the one-parameter shooting method to study numerically the problem of soliton- 



magnon scattering, as described in Ref. ^ f° r vortices in the antiferromagnet. Choosing 
the shooting parameter e, see Eq. ([To]), we "kill" the growing exponent in ( |17b|) ; as a result 
we have obtained a well-pronounced exponential decay for v m ; the corresponding oscillating 
solutions for u m differ from the asymptotical values ( |3~0"D only in the vicinity of the soliton 
core, see Fig. ||. The scattering amplitude was found from these data by comparison with 
the asymptotics (^). The results are the following: 

For all modes the scattering amplitude u m (k) tends to zero as k — > 0. In the long- 
wavelength limit the maximal scattering is related to the mode with m = —1, for which 
the behaviour of da/dk looks singular. However, other properties depend strongly on the 
soliton radius. 

In the case of large R, the scattering amplitude cr m (k) is positive for all modes in the 
long-wavelength limit. There is a maximum of <J m (k) at about kR ~ 1 for modes with 
positive m; but there are two poles in the scattering amplitude for modes with negative m, 
and one pole for m — 0. Naturally, there is no real divergence at such a pole: the physically 
observed phase of the scattering, or the phase shift 5 m = — tan cr m , varies monotonically, see 
Fig. |6|. Thus, there are principally different pictures for positive and negative m: the total 
phase shift is 



riot 




S(oo) - 6(0) = < -7T, m = , R > A . (31a) 



The scattering picture is quite different for the case of small radius solitons. Each mode 
with \m\ > 1 has a single pole only, the position k p R of which increases when m grows. For 
the special case of the mode with m — +1, the scattering amplitude is anomalously small; 
for the mode with m = — 1 there appears a second pole with anomalously small values of 
k p . These results are presented in Fig. |7|. The total phase shift is 

+7r, m > 1 

C = { °' m = 1 , , R « Ao- (31b) 

m 1 — 27r, m = — 1 

-7r, m = 0, m < —1 
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There appears discontinuous change of the total phase shift from fl31a|) to ( |31b|) for all 



modes with m ^ 0,-1 for the cases of intermediate soliton radii. This transition occurs 
independently for each mode at some critical values k c m of the wave-number. 

VI. THE SCATTERING PROBLEM IN THE FRAMEWORK OF THE 
PSEUDOPOTENTIAL METHOD 

For the analytical description of the scattering problem, let us analyze the EVP quali- 
tatively, replacing the real "potentials" fll~5|) in the EVP fll4]) by simplified ones, which are 
amenable to an exact treatment. Such an approach is well-known in the quantum physics 
of solids as the pseudopotential method. Sometimes it is useful even when the potential is 
not small, see Ref. |37| for details. 

Let us discuss the choice of the pseudopotential. The simplest way is to describe the 
soliton structure in the main approximation, using the ansatz 

r - R 

costfoO) « sign f, f = — — , (32) 

which is a good approximation for both small and large distances. Replacing the soliton 
by this configuration, which does not have any in-plane structure (sin#o — 0), reduces the 
complex EVP ( |14"D to a very simple "centrifugal model" with v m = 0: 

vy2 , . 2 . I^ + Sign £|2\ 

-V r + k H — 1 u m = 0. (33) 

This model describes quasi-free magnons in each of the regions r < R and r > R. The only 
effect of the soliton-magnon interaction is a shift of the mode indexes in comparison with 
the free magnons fl28|) : 

u (r)<xl J I™-h( At )> r<i? ' f 34 ) 

The usual matching condition for these solutions has the form 

= 0, (35) 



v! 



u 



R 



where [. . ] R = (. . -)\ R+0 — (• • -)lii-o- Calculations lead to the scattering amplitude in the 
form: 

, v = J |m-l|( X ) - J\m+l\(x) ~ J'\ m +1\( H ) [ J\m-l\(x) 

where x = kR, see Fig. |^. The total phase shift in this model 

5% = 7T • sign m. (37) 
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A. Small radius solitons, R <C Aq 



The simple "centrifugal model" with the results ( |56| ) explains qualitatively the scattering 
data for the small radius solitons, except for modes with m — 0, ±1, cf. ( |31b|) . The cause 
for these exceptions is the existence of truly local modes ("zero-ones" for m = 0, +1, and 
another one with non-zero frequency for m = —1); while such modes are absent for the 
other m. The existence of a local mode diminishes the total phase shift by 7r (and decreases 
the number of states in the continuum spectrum by 1, due to the conserved total number 
of magnon states). Such a result is well-known in Id, see e.g. Ref. |7]. Except for these 
cases, one can expect for the scattering peculiarities at/c~l/_R>>l/A . In this range the 
approximation by an isotropic magnet and BP-solitons is valid. The scattering problem for 
the isotropic magnets was studied in detail in Ref. [TP]. Over a wide range of wave-numbers, 
1 / Ao <C k <C 1 / R, the asymptotics for the Belavin-Polyakov soliton are valid,0 

rr(kR) 2 

am ^ ik) = - 2M(m + l) ' (38a) 

a -° {k) = 21n(W (38b) 
cr m= _i(A;) = n(kR) 2 In (1/kR) . (38c) 

For larger wave-numbers, k ^> Wn\/R, the scattering amplitude has a general, so-called, 
eikonal dependence, a m oc 1/kRpi 

7i(m — 1) , 

<J m {k) = 1 — ; - (38d) 



Our numerical data in this region are very close to the asymptotical behaviour fl38|) , see 
Fig. [7[ Note that the maximal scattering in the region 1/Ao <C k <C 1/R occurs for the 
rotational mode with m = 0, whose asymptotics ( |38b|) has a singularity for k — > 0. There is 
one pole only for each mode; its position corresponds to the pole in the Belavin-Polyakov 
solution, k p ~ |m|/i?S 

Let us consider the local modes, for which the total phase shift is diminished by 7r in 
comparison to the "centrifugal model" (0); it is equal to 0, — tt, —2tt for the cases m = +1, 
0,-1, respectively. 

There is one pole only for the simplest case of the mode with m = 0. The scattering 
amplitude in the range l/Ao<CA;<Cl/-Ris described by the formula (|38b| ) . The position of 
the pole is determined by the BP-result k p — 1/ The reason for such a nice agreement 
with the model of the isotropic magnet will be explained below in the text. 

There is no pole for the translational mode with m — +1. This is due to a unique 
property of the isotropic model, where such a mode does not scatter at all;El the magnon 
amplitude in this limit has the form 

which agrees with our numerical data, see Fig. This unique property of reflectionlessness 
is caused by the high internal symmetry of the problem. 
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Fig. [F](a) shows that such a property is approximately valid for the small radius soliton 
in the anisotropic magnet, too: the scattering amplitude is anomalously small, so the total 
phase shift is absent . Another effect of this symmetry in the isotropic magnet is that the scat- 
tering problem for the BP-model can be reduced to a simplified one with effectively smaller 
scattering. Such a "reduced model" has no local (and quasi-local) modes; it is described by 
shifted indexes, m — > m — lJll It is natural to use the "centrifugal" approximation for such 
a reduced model (BBJ) . As a result, a m=+ i = in the framework of the same "centrifugal 
model"; for the other modes, Eq. d37|) leads to the total phase shift 8^ = n • sign(m — 1). 
Taking into account that sign = 0, one can see that such a simple dependence for 5%? 
describes all modes, except of m = — 1. Its behaviour is most interesting. 

For the mode with m — — 1 there are two poles, while it should be one pole only 
for the Belavin-Polyakov solitons. Indeed one of the poles lies in the region where the 
Belavin-Polyakov approximation is valid; its value, k p ~ 2.79/ R ^> 1/A , coincides exactly 
with the Belavin-Polyakov dataS The second pole lies in the region k <C 1/A , where the 
approximation by the isotropic magnet is not applicable. This pole is caused by the presence 
of the local mode, the typical wave-number of such a mechanism is k p ~ >cf <C 1/Aq. The 
use of "reduced" or "centrifugal" models is certainly wrong for this case. We will discuss 
the scattering amplitude in this region later, on the basis of a more general pseudopotential 
model, see Eq. (|46a|). 



B. Large radius solitons, R Ao. Refinement of the pseudopotential model 

Let us consider the case of large radius solitons, R 3> A . Here exist local modes for 
each m. Therefore the simple "centrifugal model" as well as other simple models without 
a discrete spectrum seem to be wrong. This can be explained by the fact that in the very 
simple model (|33|) we omit the "coupling potential" ; the real potential W tends to zero at 
the origin and at infinity, but has maxima at some finite distances. 

Thus we need a more general form of the pseudopotentials which accounts for the poten- 
tials Uq, V, and W. To describe the soliton shape, we can again use the expression (|32|) for 
the out-of-plane soliton configuration. However, we need to consider the in-plane signature, 
too. According to (|32]) , it could be described as 

sm 2 e (r)— ,5(0- 

Therefore we can suppose that the "coupling potential" W(r) is nonzero only at some 
distance around the soliton radius and approximately replace it by a 5-potential as follows, 



The same assumptions lead to simple expressions for the other potentials in the EVP 



W(r) — > ^(0- (40a) 



(0), 



9 

m 



U (r) + — ±V(r 



1 - (Q/uj ) - sign g \m ± sign£| 2 c 2 

A2 + ~2 + A2^' ( 40b ) 



15 



Here c\ and c 2 are trial parameters. 

We choose this pseudopotential model (0) out of the following reasons: (i) it allows an 
exact analytical solution; (ii) it guarantees the correct asymptotic behaviour of the solutions; 
(hi) it offers a possibility to identify the parameters of the pseudopotential from the discrete 
spectrum (p3|). 

Let us consider the discrete spectrum in the framework of the pseudopotential (fiO|). 
Comparison with the eigenfrequencies (|2B[) provides a possibility to determine all trial pa- 
rameters Cfc. For the discrete part of the spectrum, the EVP (|T4|) with potentials ( f4~0| ) has 
solutions like 



u m r oc 



r oc 



Ku 



-l|l X 2 r J> 

t +ii( x r r ) 



'|m+l| 



Iro-llW r J 



r < R 

r > R 

r < R 

r > R 



x. 



^(^0 I 0±w) /-D, x± = y/(u + Q±u) /D. 



Now it is possible to calculate the eigenspectrum by matching the solutions 
( [41 b|) at the distance R, 



u 



R 



R 



C 2 

A 

c 2 

A 



A 
A 



(41a) 
(41b) 

and 

(42a) 
(42b) 



Calculations lead to the eigenfrequencies in the form (|20|) . Using the condition for the "zero- 
modes" (pl~D, the eigenfrequencies can be represented in the form (|22|); the constants ci, C2 
can be identified by comparison with the spectrum (1221), 



Cl 



P- 



a 



c 2 



— 2 + ci + ma + m 2 /5 — 



A 2 
AR 2 ' 



(43) 



The parameters a and /3 were calculated by the variational approach, a oc 1/-R 3 , /3 oc 1/R 2 , 
see Eq. ( |23a| ). 



Let us discuss the scattering problem. Solving the EVP ([[4]) with pseudopotentials (pEDP, 
we obtain for the continuum spectrum, u > ^ ga p, the following solutions. The function u 
for r < R has the form 



u(r) oc 



-f|m-i|(><2 r), when k < ^2Q/D, 
J\m-i\(k2r), when k > y/2Q/D, 



(44) 



where k^ = \/k 2 — 2Vl/ D. At large distances, r > R, it has the usual oscillating form 
(j30l). The function u has the localized form ( |41b|) . One can find the scattering amplitude 
by matching the solutions (|3H|), (0), and ( [41 b| ) at the distance R, using the condition 



[gm(fe) + ffmW] J\ m +i\{kR) - kA J' lm+ll {kR) 
kA Y'(kR) - {£ m (k)+3m(k)}Y lm+ll (kRy 



(45) 



here 
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c 2 - C, 



c 2 + x^~A 



- x+A 



K| m _!|(x+it!) 



x 2 A • 



m— 1 



X, 



J|ro-l|(**2 R) 

J\m-i\{k2Ry 



I\ m+ i\{>4R) 



k < y/2tl/D, 
k > ^2Vl/D. 



It is easy to analyze the expression for the scattering amplitude Q35D in some limiting 
cases. In the long-wavelength limit, k <C 1/R <C 1/A , one can use the asymptotics for the 
cylindrical functions J(z), Y(z) at z <^ 1 and, at the same time, the asymptotics for I(z), 
K(z) for z ^> 1. Simple calculations show that the scattering intensity is maximal for the 
mode with m = — 1, 



-(*) 



7T 



21n(l/fci2)' 



< 1, 



(46a) 



which agrees with the numerical data, see Fig. ||(b). Note, that the Eq. ( [46 a| ) was obtained 
using one condition only, kR — > 0. Therefore it is valid for the case of small R, too, explaining 
the existence of the second pole for the mode m = — 1 both for large and small radii. 

For the other values of m, we can restore the general dependence of the amplitude of the 
magnon scattering on the large radius soliton, R ^> A : 



(k) 



7T (kR/2f m+11 
\m + l|!(|m + 1| - 1) 



kR < 1. 



(46b) 



The scattering amplitude, according to Eq. (|46b|) , is positive for all modes in the long- 
wavelength limit. 

In the range l/i?C fc€ 1/V-RA the asymptotics for J(z), Y(z) for z ^> 1 are valid, 
which in the main approximation for Aq/R leads to a linear dependence of the scattering 
phase shift on kR, 



$m{k) 



\m + 1 1 7T 7T 



4 



- kR, K kR -C 




(46c) 



For this range of k, all curves look like a family of parallel lines with equal distance of about 
7r/2. Note that the real values of the trial parameters c\ and c 2 lead to the corrections in 
Eqs. (ID like A /R. 

Naturally the pseudopotential method cannot explain the scattering results for k > 1/A , 
where the influence of truly local modes and the differences between the real potential and 
the 5-function are dominant. Nevertheless, the method gives quite good results even for 
the case k < 1/Ao, where the linear asymptotics does not work. Corresponding curves are 
plotted in Fig. |S| using the general dependence fl4"5| ) with trial parameters C\ and c 2 . Fitting 
of these parameters gives a result in the zeroth-order approach for A /R. c\ ~ 0, c 2 ~ 0.4 
instead of C\ ~ 0, c 2 ~ —2, according to (H3). The difference is caused by the limitations 
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of the model fl40|). For example, the above mentioned model leads to an increase of the 
scattering phase when k increases; the total phase shift for this model is 5^ = sign m ■ it, 
as for the "centrifugal model", see (|3"T|), i.e. the qualitative properties of these models are 
equivalent. The appearance of a local mode diminishes the total phase shift by 7r. As a 
result, there is no total phase shift for all modes with positive m; but 5 tot = —2n for all 
modes with m < 0. 

The scattering picture in the short-wavelength limit, k ^> \m\/R, gives the general, so- 
called eikonal dependence o m ~ A m /k. However, in contrast to the case of the isotropic 
magnet ( |38d| ), the scattering amplitude has the same sign A m < for all modes due to 
the presence of the effective attractive potential, which causes the existence of local modes. 
When the soliton radius decreases and the local mode disappears, the sign of a m changes; 
the result becomes similar to 



Thus, we explain the particularities in the scattering picture of EA magnets, in com- 
parison with the isotropic case, by the influence of local modes. To examine this influence, 
let us consider the scattering data a m for the modes m = +2 and m = —2 with different 
soliton radii R as functions of the parameter kR, i.e. a = cr m (kR); this gives a possibility 
to compare the scattering data with the isotropic model, see Fig. [10]. For the small radii 
the scattering picture is similar to the BP-case (solid lines in Fig. ITUf). When the soliton 
radius increases, the scattering phase jumps up by tt, decreasing the total phase shift. Such 
a bifurcation takes place for the mode with m = +2 at R = 2.24A , and for the mode with 
m = — 1 at R — 1.52Ao- These values agree with the crossover points R c ml see Table ||. 
Thus, it is natural to suppose that the bifurcation is connected with the appearance of the 
quasi-local modes inside the continuum spectrum for R > R c m . Simple calculations show 
that the positions k p of the poles correspond to the frequencies of quasi-local modes, which 
can be estimated by Eq. ( |23b| ). 



VII. TRANSLATION AL MODES AND DESCRIPTION OF THE SOLITON 

DYNAMICS 

The foregoing analysis of the spectrum for the modes with |m| = 1 can be used to describe 
the translational motion of the soliton. Such an approach was first proposed in Ref. |38|, where 
an effective mass of the vortex in a 2d EP magnet was calculated from numerical data for 
the eigenmode spectrum of the finite size magnet. An important progress for the vortices in 



FM was achieved in Ref. ETt the specific structure of the spectrum of modes with \m\ — 1 



(namely the existence of a Goldstone mode with u oc 1/L 2 , and a number of doublets 
with mean frequency To oc 1/L and a small splitting) leads to non-Newtonian equations of 
motion with 3 rd , 5 th , ... time derivatives. Such equations were derived and discussed earlier 



in Ref. £9],using a phenomenological approach, for a review see Ref. [E]. The analytical 
calculationEJ verified their adequacy; besides, it allows to calculate the coefficients in the 
equations of motion with an error of only about 0.8%. Thus, it is possible to describe in this 
way such unusual properties as the appearance of the non-Newtonian equations of motion. 

The spectrum of magnon modes with \m\ — 1 has a different structure for the isotropic 
FM. First of all, there are no doublets]^ for an approximate description of the dynamics only 
two modes are sufficient: the translational Goldstone mode with dependence u oc 1/L 2 and 
another mode with m = — 1. Therefore, the effective equations of motion have a Newtonian 
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form with an effective mass, which diverges as L 2 , see Ref. |19|. Such behaviour agrees with 
the direct calculations done by Zaspel.ta 

Let us note that all magnets mentioned above, i.e. EP and isotropic systems, have a 
gapless dispersion law and, as a result, a strong interaction of the soliton with the boundaries 
of the magnet; the interaction force is proportional to 1/L for the EP vortices^ and 1/L 2 for 
the Belavin-Polyakov solitons.0 The absence of a gap manifests itself in the particularities 
of the dynamical coefficients (the effective mass M for the BP-soliton in the isotropic FM0 
and the coefficient G 3 in the term with [e z x d 3 X/dt 3 ] for the vortex in the EP FmH). 
This can be explained as follows. The existence of the mass M and the gyroscopical force 
G [e z x dX / dt] leads to a finite frequency of the Larmor precession of the soliton. For a 
gapless dispersion law, this frequency lies in the continuous spectrum. Thus, the Larmor 
motion of the vortex leads to the generation of magnons. For a finite-size system without 
dissipation, such magnons are distributed through the whole magnet; the Larmor dynamics 
of the soliton centre is strongly coupled with a "magnon cloud", having a scale of about 
the system size. Therefore, the position of the soliton centre, which can be determined as a 
point with m z = — 1 (BP-soliton), or m z = ±1 (vortex), in fact plays the role of a collective 
variable which governs the motion of the "magnon cloud" . Then it is not surprising that the 
corresponding equations of motion are nonlocal, leading to the divergence of the dynamic 
coefficients as L —>■ 00. 

The situation is quite different for the case of the EA FM. Our numerical analysis and 
analytical calculations show that the picture of doublets for local modes with \m\ and — \m\ 



is valid for large R, see Sec. [TV]. On the other hand, the eigenfrequencies of these modes lie 
in the gap for the large soliton radii, u < u gauP , so one can expect the existence of a finite 
mass in the effective equations of the soliton motion. Therefore for the phenomenological 
description of the soliton dynamics we can use the 2 nd -order differential equation 



dX 



F e (X) , (47) 



which corresponds to the picture with the lowest doublet of local modes. Note that further 
generalizations for the finite size magnet take into account the next doublets of quasi- 
continuous modes. Here the hierarchy of the effective equations of motion containing only 
even-order time derivatives manifests itself. In Eq. (f4~7|) X describes the position of the 
soliton, M is the mass coefficient, G the gyrocoefncient, e z the unit vector along the easy 
(z) axis, and F e the external force acting on the soliton due to the boundary or other 
solitons. Assuming that deviations of the soliton from the equilibrium position are small, 
F e « —aX, the effective equation of motion ( fT7| ) can be solved by a harmonic ansatz, which 
leads to 



G G 2 a 

^ 1 = 2M ± V4i^ + M- (48) 

For the localized solitons, the interaction between soliton and boundary (or between two 
solitons) has an exponential decay as exp(— L/r ), r = a/ D/uj ga , p , see Eq. (|9|). Thus ot/M 
can be neglected for a large system size L, so the frequencies ([Hf) have the form: 

u = 0, wi = G/M. (49) 
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The zero frequency uo corresponds to the Goldstone mode (position shift of the soliton), in 
which all spins rotate with the frequency of the soliton precession Q in the laboratory frame. 
The presence of the rotating frame frequencies in the phenomenological description of the 
soliton dynamics has a simple explanation: 

The frequencies ( pE9|) describe, in essence, the motion of the soliton centre. Let us discuss 
their link with the frequencies of small oscillations, see (|T3|) . It is easy to write down 
expressions for the magnetization components in the laboratory coordinate frames: 



m x + im y = j^sin 9 + cos 9 (u a + v a ) cos(mx + uj a t) 

a 

(u a -v a )sm(mx + 0Jj)] • e iq * +ia \ 
cos6> — sin# {u a + v a ) cos(mx + Cj a t). 



m 2 



One can see that the in-plane components have no simple frequency dependence, including 
the combined values u ± IX At the same time, the out-of-plane components, m z , depend 
on the frequency Co in the rotating frame only. In our phenomenological description only the 
^-components contain information about the position of the soliton centre X, determined 
by the condition m z = —1. Thus, just the frequency u, which is calculated above and 
plotted in all figures, determines the soliton motion. 

The solution of the effective equation ( pET| ) with the nonzero frequency uji allows us to 
calculate the effective mass of the soliton, M = GfiO\. The value of the gyroconstant is 
well-known, G = 4nqA/D, see Refs. [£3],[7|. As a result, the soliton mass is finite for every 
soliton radius (the local mode with u% exists for every R, see Fig. |]). This fact corresponds 



qualitatively to the calculations in Ref. §1]. The soliton mass we have obtained numerically 



from the c^i-data is plotted in Fig. O as a function of R. 



In the limit of a large radius soliton the formula ( p3d|) yields the dependence 

M = M - , M = -^, i?»A , (50a) 



which agrees with the results obtained in Ref. [41] and our numerical data, see Fig. |12. 
In the case of small radius, R <^ A , the dependence (|26|) is valid, thus 



M = M .^- = ^-rl J R«A , (50b) 

where r was introduced in Eq. (|^). This result corresponds to the one for the Belavin- 
Polyakov soliton,0 M = 4iiAD~ 2 ■ L 2 , because here the characteristic length r is to be 
replaced by the system size L. 



VIII. CONCLUSION 

We have studied the magnon-soliton system in the model of the 2d Heisenberg easy- 
axis ferromagnet. Combining numerical and analytical methods, we have obtained complete 
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results about bound and scattering magnon states in the system. A rich spectrum of truly 
local modes was found along with bifurcations of these modes with the change of the soliton 
radius. For the modes with higher azimuthal numbers m, we have verified the picture 
of doublets with small splitting for large soliton radii R. The spectrum changes strongly 
with decreasing of R: all modes with \m\ > 1 "leave" the region of the discrete spectrum, 
transforming to quasi-local modes. Such modes could be observed experimentally by soliton 
magnetic resonance as was done for Id solitons.i 

The influence of truly local modes is important for the scattering problem, leading to 
the bifurcations of the phase shift. As a result, when R < 1.52A , the scattering picture is 
qualitatively the same as for the isotropic magnet which could be explained in the framework 
of a simple "centrifugal model". There is one exception only, m — — 1, where the local mode 
exists for every soliton radius. 

Our investigations can be applied to the description of the eigenmodes of magnetic dots. 
In the present paper we develop the theory of local modes on the soliton background; it 
could be a good guide for the study of the normal modes in the vortex-state magnetic dots. 
Our theory is constructed for the soliton in the model (|1|), where the soliton is stabilized by 
the internal precession. It is clear that such a model cannot guarantee the quantitative cor- 
respondence with the case of vortex-state magnetic dots, where the static soliton structure 
is stabilized by the magnetic-dipole interaction. We did not consider this type of interaction 
in the paper, as it is difficult to account for. Nevertheless we believe that the main features 
of the problem studied above are generic. First of all, we expect the appearance of modes 
with anomalously small frequencies, e.g. the mode of the translational oscillations of the 
vortex centre. The non-zero frequency of this mode is caused by the interaction with the 
boundary only. Secondly, doublets with m = ±n, n > 1, should appear, when the radius of 
the vortex is rather large. 

We have used the results on the local modes for the most interesting case of the trans- 
lational modes with \m\ = 1 to describe the soliton motion in the infinite-size magnet: it is 
possible to identify the soliton mass, which is finite due to the localized soliton structure. In 
contrast to both the Belavin-Polyakov soliton in the isotropic FM and the magnetic vortex 
in the easy-plane FM, the soliton motion in our case of the easy-axis FM is similar to the 
motion of a finite mass charged particle in a magnetic field. 

Thus, to the best of our knowledge, the localized precessional soliton in the easy-axis FM 
is the first example of a 2d topological magnetic soliton with truly particle-like properties. 
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TABLES 



TABLE I. Fitting of the trial solution (|l0|) to the exact numerical solution for 9o(r). Radius of 
the soliton is defined from the condition Oq(R) = tt/2. 



Frequency of 


Soliton 


Fitting para- 


Std. error. 


precession, Q/u>o 


radius, R/Aq 


meter, i?( fit ) / A 


xl0~ 2 


0.05 


20.38 


20.38 


0.6 


0.10 


9.95 


9.94 


1.2 


0.15 


6.58 


6.58 


1.9 


0.20 


4.88 


4.87 


2.6 


0.25 


3.82 


3.81 


3.2 


0.30 


3.09 


3.07 


3.6 


0.40 


2.07 


2.07 


4.0 


0.50 


1.35 


1.37 


3.8 


0.60 


0.81 


0.85 


3.2 


0.70 


0.42 


0.44 


2.2 


0.80 


0.15 


0.16 


1.1 


0.85 


0.07 


0.07 


0.9 



TABLE II. Critical parameters for some lower modes with \m\ > 1. Remember that local 



modes always exist for m = 0, ±1. 



rn 










data 


data 


from Eq. flU) 


-2 


1.52 


0.47 


0.45 


+2 


2.24 


0.38 


0.38 


-3 


2.73 


0.33 


0.30 


+3 


3.22 


0.29 


0.28 
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FIGURES 




FIG. 1. Shape of 6 (r) for different soliton radii R/A (x = 0.4, ■ = 5, = 10). Lines: 
numerical data, symbols: fitting data with the trial function (|i~0|) 




0.5 1 1.5 2 2.5 3 3.5 4 
Distance, r/R 

FIG. 2. Shape of the "coupling" potential for different soliton radii 
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(a) 

1 | , . . r 




Soliton radius, R/A, 



(b) 




Frequency of soliton, Q./co 

FIG. 3. Eigenfrequencies of local modes as functions of: (a) soliton radius and (b) frequency 
of internal precession of the soliton. Lines: theoretical results from Eqs. ( |23b| ) and (|27|); symbols: 
numerical data. 



26 



(a) R=10A, 



1.2 




1 2 3 4 5 
Distance, r/A 

FIG. 4. Magnon amplitudes for local modes with m = — 1 
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Soliton shape 
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FIG. 5. Oscillating magnon amplitudes u m for soliton with radius R = 10A for different wave 
numbers k 
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(a) m>0 
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Wave number, kA 
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(b) m<0 
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Wave number, kA 

FIG. 7. Scattering for a small radius soliton, R = 0.2Aq. Solid line: analytical asymptotics 
for k >C 1/Ao for the mode m = — 1 from Eqs. ( 46a ); dashed lines: analytical asymptotics for 
k 3> 1/Aq from the Eqs. (p3); symbols: numerical data 
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FIG. 8. Scattering for the "centrifugal model' 
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(a) m>0 




012345678 
kR 

FIG. 9. Scattering amplitude S m for large radius soliton (R = 10Ao) vs. kR in the range 
kAo < 0.8. Lines: analytical solution from (|45|); symbols: numerical data 
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FIG. 10. Bifurcations for the scattering data 8 m vs. kR: (a) critical soliton radius 
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1.52Aq. Solid line corresponds to the Belavin-Polyakov asymp- 
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1 2 3 4 5 6 7 
Distance, r/A 

FIG. 11. Magnon amplitude u for the translational mode with m = +1 (R = 0.2Ao, fcAo = 3). 
Lines: exact analytical solution (39) for the isotropic magnet; symbols: numerical data for the EA 
FM 
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FIG. 12. The R dependence of the mass. Lines: analytical asymptotics from Eqs. (|50|); symbols: 
numerical data 
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